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ABSTRACT 

The  use  of  an  electronic  analog  computer  in  solving  the  differential 
equation  of  motion  of  the  free  piston  gas  generator  has  been  investigated. 
To  accomplish  this  end,   the  free  piston  has  been  treated  as  a  specialized 
mass-spring  system,  wherein  the  forces  were  derived  from  thermodynamic 
considerations  rather  then  from  elastic  considerations.      The   first  law  of 
thermodynamics  was  solved  on  the  computer  three  concurrent  times  in  order 
to  produce  electronic  analogs  of  an  internal  combustion  engine,   a  recipro- 
cating compressor,   and  a  gas  spring.     Coulomb  and  viscous  friction  were 
also  simulated.      The  results  of  these   four  analogs  were  then  combined  and 
integrated  as  indicated  by  Newton's  second  law.      The  computer  yielded  a 
solution  which  represented  the  piston  motion0     A  total  of  26  amplifiers, 
seven  diode  limiters,   and  six  electronic  multipliers  were  used  to  solve 
the  problem o 

The  characteristics  of  an  existing  free  piston  machine  were  intro- 
duced into  the  computer  as  a  final  check  of  the  analog.      The  accuracy  of 
some  of  the  results  was  less  than  desired  due  to  errors  in  a  two-quadrant 
electronic  multiplication  involved  in  the  gas  spring  analog.      The  other 
three  analogs  performed  very  satisfactorily.     For  accurate  and  repro- 
ducible results,   some  other  method,   such  as  direct  function  generation 
or  single  quadrant  multiplication,  must  be  used  in  the  simulation  of 
the  gas  spring. 
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1.   Introduction, 

A  groat  deal  of  offort  has  been  put  into  the  solution  of  the  free 
piston  engine  problem.  All  of  the  methods  currently  in  use  are  digital 
methods  involving  a  point-by-point  calculation  of  the  variables  asso- 
ciated with  such  a  deviceo   These  methods  are  time-consuming  and  cumber- 
some due  to  both  the  trial  and  error  element  and  the  many  geometric  and 
operating  parameters  involved.   In  addition,  the  digital  methods  can  be 
made  tractable  only  by  incorporating  many  numerioal  assumptions j  some 
of  which  introduce  considerable  error. 

London    has  developed  a  " Thermodynamio-dynamic  analysis"  whereby 
the  calculated  parameters  of  a  "standard"  machine  may  be  scaled  by 
certain  affinity  relationships  to  establish  the  design  of  larger  or 
smaller  machines.  This  method  eases  somewhat  the  numerioal  burden  of 
calculation,  but  the  results  proved  to  be  quite  sensitive  to  the  assump- 
tions madej  and  in  order  to  predict  the  characteristics  of  an  existing 
machine,  the  assumptions  were  relaxed  in  four  successive  groups  before 
the  desired  aocuraoy  was  obtained.   This  procedure  still  presented  a 
formidable  numerical  task. 

In  addition,  a  solution  by  digital  methods  becomes  more  oomplex  in 
the  investigation  of  part-load  characteristics. 

The  purpose  of  this  investigation  is  to  determine  the  possibility 
of  solving  the  differential  equation  of  motion  of  the  free  piston  on 
the  Boeing  Electronic  Analog  Computer  (BEAC),  whioh  is  described  in 
Appendix  VI,  If  this  can  be  done,  the  computer  will  electronically 
simulate  the  free  piston,  presenting  as  a  solution  the  piston  motion, 
with  most  of  the  pertinent  quantities  also  available.  The  advantages 
of  this  method  inoludei 

1,   The  feedback  prinoiple,  typical  of  analog  computation, 
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eliminates   trial  and  error. 

2,  All  of  the  parameters  may  be   taken  into  consideration,   and 
the  effeot  of  each  may  be   observed, 

3,  The   solution  is   immediate,   thus  permitting  the   study  of  many 
combinations  of  parameters  in  a  short  time, 

4,  Fewer  simplifying  assumptions  are  required,  and  these  may  be 
relaxed  to  any  degree   compatible  with  the  size  of  tiie  computer, 

5,  Part-load  behavior  can  be   studied  as  easily  as  full-load 
behavior, 

6,  Time  may  be   scaled  to  permit   slow-motion   studies. 


2.  General  Background  of  Free  Piston  Systems. 

Figure  1  illustrates  schematically  two  of  the  possible  free  piston 
configurations.  Figure  1(a)  illustrates  the  free  piston  air  compressor, 
where  the  useful  output  is  compressed  air)  and  Fig.  1(b)  illustrates  the 
free  piston  gas  generator,  where  the  useful  output  is  high  pressure- 
temperature  gas  for  use  in  a  turbine. 

This  investigation  will  encompass  only  the  free  piston  gas  generator. 
A  slightly  larger  computer  is  required  to  simulate  the  free  piston  com- 
pressor. 

The  free  piston  machine  depicted  in  Fig.  1(b)  consists  of  two 
opposed  pistons,  the  inboard  ends  of  whioh  operate  as  a  two  stroke  diesel 
engine.  Another  part  of  each  piston  acts  as  a  reciprocating  air  com- 
pressor. As  the  pistons  move  outward,  the  expansion  work  of  the  engine 
and  compressor  is  stored  as  energy  in  the  bounce  chambers.   This  energy 
is  then  used  on  the  inward  stroke  for  the  compression  and  delivery  of 
the  air  in  the  compressor  and  the  compression  of  the  engine  oharge. 

The  outstanding  features  of  this  free  piston  machine  aret 

1,  Simplicity  of  construction  due  to  absence  of  crank 
meohanisms, 

2.  Opposed  piston  configuration  eliminates  vibration. 

3#  No  piston  side  thrust  arising  from  crank  inertia  forces. 
In  the  case  of  the  gas  generator,  it  would  be  well  to  review  the 
features  of  a  typical  cycle. 

Referring  to  Fig.  1(b),  which  illustrates,  say,  the  pistons  just 
starting  the  instroke.  Beginning  at  this  point,  the  following  successive 
processes  may  be  visualized: 

1.  High  pressure  in  the  bounce  space  forcing  the  pistons 
together;  oompressing  and  delivering  the  air  in  the  compressor  spaoe,  and 
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Figure  1.  Typical  Free  Piston  Configurations. 


oompressing  the  engine  oharge  when  the  scavenge  ports  are  oovered. 

2.  Reversal  of  piston  motion  occurring  during  injeotion  and 
burning  of  fuel,, 

3.  The  expansion  of  the  fluid  in  the  engine  cylinder  and  the 
intake  stroke  of  the  compressor;  while  the  bounce  cylinder  is  storing  the 
expansion  energy  of  the  engine  and  compressor, 

4.  Opening  of  the  engine  ports  during  the  latter  portion  of 
the  outward  stroke,  delivering  working  fluid  to  turbine. 

5*  Pressure  in  bounce  space  causes  reversal  of  piston  motion, 
and  cycle  repeats. 

The  preceding  cycle  is  represented  on  Temperature -Entropy  coordi- 
nates in  Fig.  2.   It  is  to  be  noted  that  the  processes  involve  different 
masses,  and  can  be  represented  only  on  a  "per  pound"  basis.  This  diagram 
is  similar  to  that  in  common  use  in  free  piston  literature.  It  should  be 
pointed  out  that  the  following  two  assumptions  are  made  in  producing  this 
diagrami 

1.  The  scavenge  air  is  furnished  at  oonstant  pressure  and 
temperature,  related  to  the  compressor  d is charge  state, 

2.  Delivery  to  the  turbine  occurs  at  constant  pressure  and 
temperature. 

These  two  assumptions  are  intimately  connected  with  the  mechanical  design 
of  the  engine  ports  and  the  physical  size  of  the  turbine  receiver.  It 
is  beyond  the  scope  of  this  investigation  to  consider  suoh  problems.   The 
foregoing,  therefore,  constitute  the  first  two  assumptions  in  this 
investigation.   It  must  be  pointed  out,  however,  that  these  assumptions 
must  be  relaxed  if  it  is  desired  to  investigate  transient  behavior,  inas- 
much as  these  particular  parameters  are  the  primary  transients.   Before 
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Figure  20     Temperature  -  Entropy  Representation  of  Fret  Pist  n 
Generator-Turbine  System0 


proceeding  further,  the  following  are  some  additional  assumptions 
covering  the  device  as  a  whole: 

1.  That  the  compressor  intake  is  at  70°F.  and  14.7  psia. 

2.  That  the  operating  fluid  obeys  the  ideal  gas  laws* 

3,  A  nominal  pressure  drop  of  five  per  cent  across  each  valve 
and  engine  port.   This  is  arbitrary,  and  may  be  refined  to  any  desired 
value  • 

4,  That  the  geometry  is  symmetrical  with  respect  to  a  plane 
passing  through  the  center  of  the  machine  perpendicular  to  the  hori- 
zontal axis.   The  end  projection  of  this  plan©  will  be  referred  to  as 
the  "center  line." 


♦Refer  to  Appendix  I  for  all  of  the  assumptions  used, 


3»   The  General  Analog  of  the  Free  Piston. 

From  a  dynamic  point  of  view,  the  piston  in  a  free  piston  machine  may 
be  considered  as  a  rather  special  sort  of  mass-spring  arrangement,  wherein 
the  "springs"  yield  forces  derived  from  thermodynamic  considerations 
rather  than  from  elastic  considerations. 


Figure  3.  Spring  Analogy  of  the  Free  Piston. 


The  system  may  be  visualized  as  similar  to  that  depicted  in  Fig.  3  above, 
where  the  spring  forces  are  proportional  to  the  conventional  Pressure- 
Volume  representation  of  the  appropriate  space.  Not  shown,  but  also 
present,  is  the  foroe  of  static  and  viscous  friction.  Qualitative 
representations  of  these  forces  appear  in  Fig,  4. 

To  illustrate  the  analog  computer  approach  to  this  tyoe  of  problem, 
consider  the  simple  mass-spring  arrangement  illustrated  in  Fig.  5.   The 
differential  equation  of  this  motion  may  be  immediately  written  down  as 
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Figure  5.   Simple   Spring-Mass   System, 


where  m  is  the  mass 

v  is  the  damping  oonstant 

k  is  the   spring  oonstant 

^  is  the  acceleration 

7  is  the  velocity 

f  is  the   displacement 


By  solving  for   the  highest 
derivative, 


f-W£? 


V  /7n. 


Recalling   that  the   operational  amplifier   oan  add  and/or   integrate 
with  respect  to  time,   a  computer  soheraatic  may  be   immediately  devised  to 
solve   this  problem.      This  is  illustrated  below  in  Pig,  6. 
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Figure  6,  Electronic  Analog  of  the  Simple  Spring-Mass  System, 


•Refer  to  Appendix  II  for  analog  computer  table  of  symbols. 
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This   simple  problem  points  out  the   utility  of  analog  feedbaok.      That  is 
to  say,    the   immediate   stabilization  to   the  correct  solution  by  using  this 
solution  to  generate  the  derivatives,  which  are  then  integrated  to  form 
the   solution  again. 

Turning  now  to   the  free  piston  problem,   and  referring  to  Pig,   3;    the 
equation  for   this  motion  may  be  formulated  in  accordance  with  Newton's 
second  law  as 

MX    *  **(*,%)   +    FC(X,X)   ~  FbU,i)   "  FF(  i   )•  3ol 

where  U   is  the  piston  mass 

Fg'z/L)  is  the  force  on  the  engine  piston 

F'o(z,^)  is  the  force  on  the  compressor  piston 

F,(^X)  is  the  force  on  the  bounce  piston 

Fj^)  is  the  force  due  to  friction 

x  is  the  piston  displacement 

%  is  the  piston  velocity 

X  is  the  piston  acceleration 

In  a  manner  similar  to  the  simple  mass- spring  problem,  a  computing  ar- 
rangement can  be  immediately  proposed,  Suoh  an  arrangement  is  outlined 
in  Fig,  7.  In  this  figure  the  blook  representing  the  engine,  compressor, 
bounoe,  and  friotion,  are  rather  complex  computing  networks  with  many 
more  inputs  than  shown.   The  overall  principle,  however,  is  basically 
the  same  as  that  outlined  in  Fig,  6, 

This  fact  has  two  very  important  consequences  in  that  : 

1,  A  completely  compatible  set  of  parameters  must  exist  in 
order  for  the  computer  to  produce  a  sustained  oscillation.  If  this  con- 
dition is  not  met,  a  situation  exists  which  is  analagous  to  overdamping 
or  negative  damping  of  the  simple  system, 

2,  Sinoe  equation  3,1  also  expresses  the  work-energy  equality, 
a  solution  of  this  equation  will  automatically  satisfy  the  condition 
that  the  engine  work  appears  as  the  sum  of  compressor  and  friotion  worko 

The  computer  arrangement  as  shown  in  Fig,  7  will  produce  displace- 
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ment  as  a  function  of  time.  Also  available  will  be  the  various  forces 
or  pressures  contributing  to  the  motion.  Related  items  suoh  as  air-fuel 
ratio,  scavenge  ratio,  and  power  output  will  not  be  available  directly 
from  the  oomputer  but  oan  be  easily  calculated. 
The  inputs  to  the  system  will  be  : 

1,  Geometric  parameters,  which  will  be  fixed  for  a  given 
maohine. 

2.  Initial  pressure  conditions,  which  aret 

a.  Intake  and  delivery  pressure  of  the  compressor, 

b.  Pressure  at  the  start  of  engine  compression, 
o«  An  initial  pressure  in  the  bounce  space, 

3,  Amount  and  heating  value  of  fuel. 

4.  Thermodynamic  parameters  suoh  as  poly tropic  exponents, 
and  ratio  of  specif io  heats. 

All  of  the  above  parameters  may  be  varied  to  suit  almost  any  con- 
dition compatible  with  the  assumed  geometry.  For  instanoe,  part  load 
behavior  may  be  simulated  by  merely  changing  the  appropriate  pressure 
initial  conditions,  and  then  determining  the  range  of  fuel  energies  for 
which  the  solution  is  stable.  Of  course,  a  different  amount  of  fuel 
applied  to  the  same  pressure  oonditions  will  lead  to  a  different  dis- 
placement-frequency solution. 

The  usefulness  of  the  analog  computer  approach  to  this  problem  is 
thatt 

1.  The  parameters  are  introduced  as  various  constants  and 
initial  conditions  and  are  not  manipulated  at  all. 

2.  The  oomputer  solves  the  differential  equation  of  motion. 

3.  Other  quantities  may  be  determined  by  a  minimum  amount 
of  calculation, 

13 


It  is  to  bo  re-emphasized  that  the  analog  computer  eliminates  the 
trial-and-error  elements  from  this  problem  if   the  geometry  is  used  as  a 
starting  point. 
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4.  Geometrio  Parameter 8, 

There  are  eight  primary  geometric  parameters  and  variables  asso- 
ciated with  the  analog  approach  to  the  free  piston  problem*   These  aret 

1,  Area  of  engine  piston,  Ag, 

2,  Area  of  compressor  piston,  AQ. 

3,  Area  of  bounce  piston,  A^. 

4,  Engine  piston  displacement  from  oenterline,  x» 

5,  Compressor  piston  displacement  from  head,  y. 
6*   Bounoe  piston  displacement  from  head,  z, 

7.  Engine  piston  position  at  port  closure,  x«# 

8.  Piston  mass,  M. 

Referring  to  Fig,  1(b),  it  may  be  readily  seen  that 

Ab  "  As  "h  Ac*  4-x 

y  ■  x  +  o,  4.2 

and 

z  s  b  -  Xj  4,3 

where  o  and  b  are  the  respective  clearances  of  the  compressor  and  bounoe 

oylinders  when  the  engine  pistons  are  touohing.  Since  -the  piston  is  in 

one  piece,  it  is  obvious  that 

dx,  ..  dy  _    dz  4,4 

dt  ~  "3$  ~  "  "dT  *  ' 

These  relationships  are  necessary  in  that  x,  y,  and  z  and  their  time 
derivatives  are  used  in  the  generation  of  the  appropriate  forces. 

The  mass,  M,  may  be  taken  as  any  reasonable  value  compatible  with 
piston  size,  A  computer  solution  indicating  a  piston  density  greater 
than  that  of  steel  must  be  rejected. 


16 


The   treatment  of  the  free  piston  generator  will  encompass  one-half 
of  the  system,  which  will  be  the  right  hand  piston  as   Illustrated  in 
Fig.  3.     Positive  forces  will  be  considered  to  be  aoting  to   the  right. 
Since  one  of  the  computer  processes  will  be  the   summation  of  these  forces, 
the  associated  sign  inversion  requires   that  the  negatives  of  the  engine 
and  compressor  forces  be  generated.     Since  one-half  of  the  system  is 
being  treated,  all  extensive  quantities  must  be  multiplied  by  two. 
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5o  Thermodynamic  Parameters. 

Inasmuch  as  the  computer  is  solving  a  mechanical  problem,  all  thermo- 
dynamic quantities  will  be  assumed  to  be  given  in  compatable  mechanical 
units.  The  use  of  Joule's  constant  will  always  be  implied  whenever 
necessary.  Energy  leaving  a  system  as  heat  will  be  considered  negative. 
5<>1  Pressures. 

The  pressure  taken  as  the  fundamental  variable  will  be  the  turbine 
receiver  pressure,  P^.  This  will  be  assumed  to  remain  cyclically  constant. 
Using  this  as  a  starting  point,  it  is  possible  to  work  backwards  through 
the  system  to  establish  all  of  the  other  pressures.  These  are 

5.1.1 
5.1.2 
5.1.3 
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where  Pp  is  the  engine  pressure  during  scavenging, 
Ps  is  the  pressure  in  the  scavenge  cylinder, 
P<j  is  the  compressor  delivery  pressure, 
P^  is  the  compressor  intake  pressure, 
P0  is  the  ambient  pressure,  and 
r  is  the  ratio  of  downstream  to  upstream  pressure  for 

"e"  -  Engine  exhaust  ports, 

"s"  -  Engine  scavenge  ports, 

"dH  -  Compressor  discharge  valves,  and  / 

"i"  -  Compressor  intake  valves, 

5.2  Volumes. 

In  all  cylinders,  the  volume  is  the  product  of  total  displacement 
and  area. 

5.3  Temperatures.  ' 
Temperature  will  not  be  used  as  such,  but  rather  indirectly  in  the 

determination  of  the  mass  of  air  charged  to  the  engine  and  the  power 
output  of  the  machine. 
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5*4  Mass  of  Engine  Charge* 

The  mass  of  air  charged  to  the  engine  is  not  at  all  a  necessary  para« 
meter  insofar  as   the  computer   is  concerned,  but  must  be  known  in  ortfer  to 
determine  the  air-fuel  ratio. 

Using  the  assumption  that   there   is  no  heat  gained  or  lost  in  the 
scavenging   space,  it  can  be  postulated  that,  upon  completion  of  the  sca- 
venging process,    "the  engine   charge   temperature,   Tp,  will  be  essentially 
that  temperature  of  the  compressor  discharge.      The  mass  of  engine  charge, 
me,  will  then  be 

me  s  Pp  Vp/RTp,  5.4.1 

where  Vp  is   the  engine  volume  at  the  point  of  port  closure. 
Applying  the  compression  ratio  of  the  compressor  to  the  above  equation 
yields 


ML 


'  PA1**'1    '  5*4"2 

W 

where  T0  is  the  ambient  temperature 

nc  is  the  polytropio  exponent  of  compression  in  the  compressor. 

This  result  may  be  expressed  in  terms  of  P^  by  applying  relationships 

5,1.1  and  5,1.6! 

*e»   ^Jl (AjgJ  ^   C)        ,         6.4.3 

RX  re 

If  all  r's  are  assumed  the  same,  the  above  equation  reduoes  to 

VP 


Me  - 


ft?: 


-I 


5,4  04 


Assuming  that  r  is  0,95,  and  further  postulating  that  nQ  is  1,30,  the 
above  relationship  becomes 

■"%    =     0.0113  Pt°'\  «.4.» 
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whore  P^.  la  expressed  in  pounds  per   square   inch.     For  convenience, 

equation  5.4J5   is  presented  graphically  in  Fig.  8, 

5,5  Power  Output  of  the  System. 

To  determine  the  power  delivered  to   the  turbine,  an  overall  energy 

balance  is  necessaryt 

mdh0  -}-  mf  H7  -  (mf  +  md)ht  f  Qr9j,  5.5,1 

where  md  is  the  mass  of  air  delivered  by  the  compressor  per  cycle, 
mf  is  the  mass  of  fuel  injected  per  cycle. 
h0  is  the  enthalpy  of  entering  air. 
ht  is  the  enthalpy  of  exit  gases. 
BV  is  the  lower  heating  value  of  fuel 

Qrej  is  the  energy  as. heat  leaving  the  system  via  the  cooling 
water  per  cycle. 

Equation  5.5.1  makes  the  customary  assumption  that  the  sensible  enthalpy 
of  the  fuel  may  be  negleoted. 

The  quantity  md  may  be  readily  determined  by  noting  from  the  com- 
puter solution  the  length  of  the  intake  portion  of  the  compressor  stroke, 
and  applying  the  ideal  gas  law 

^d  =  piA0Yi/ RT0,  S        5.5.2 

where  Y^  is  the  length  of  the  intake  portion  of  the  compressor  stroke. 

The  quantity  mf  HV  is  a  known  operating  variable  and  1^  is  tabulated. 
The  quantity  Qrej  is  composed  of t 

1.  Energy  as  heat  leaving  the  engine. 

2.  Energy  as  heat  leaving  the  compressor. 

3.  Energy  lost  through  frictional  heating.  It  will  be  assumed 
that  this  will  be  totally  lost  to  the  coolant. 

The  determination  of  Qr9j  will  be  deferred  until  more  speoifio 
differential  equations  are  formulated,, 

Assuming  for  the  moment  that  Qj^j  is  known,  the  only  remaining  un- 
known in  equation  5.5.1  is  h-t;#  the  enthalpy  of  the  exit  gasses.  From 
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this  quantity,   the  temperature  of  the  exit  gas  may  be  determined,  and 
also   the  ideal  isentropic  work  available  per  pound  of  fluid.      The  pro- 
duct of  this  work,   the   total  mass  delivered  per   cycle,   and  the  frequency 
will  then  give  the  air  horsepower  of  the  system.      It  is   to  be  noted  again 
that  this   treatment  involves  one-half  of  the  machine,   so   that  the  "total" 
mass  mentioned  above   is  twice   that  given  in  equation  5.5.2. 
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6.  Electronic  Analog  of  an  Internal  Combustion  Engine. 

6.1  General, 

The  force  exerted  upon  the  engine  piston  may  be  most  appropriately- 
represented  by  referring  to  a  typical  indicator  diagram  such  as  depicted 
in  Fig.  4(a)  on  page  9.  Since  the  free  piston  device  has  no  crank  mech- 
anism, the  only  fixed  geometrical  point  is  that  of  port  closure.  The 
thermodynamic  state  at  this  point  will  be  known  if  it  is  assumed  that 
successive  cycles  are  identical.  All  other  point  of  force  as  a  function 
of  displacement  are  determined  by  the  dynamics  of  the  system.  In  order 
to  simulate  this  situation  it  is  necessary  to  obtain  a  computer  solution 
of  the  differential  equation  of  the  first  law  of  thermodynamics 0 

In  the  following  analysis,  the  approximate  air-standard  is  used,  which 
assumes: 

1.  Air  is  the  working  fluid,  and 

20  The  mass  of  fluid  within  the  cylinder  remains  constant. 

The  latter  assumption  is  also  derived  from  the  following  two 
practical  considerations : 

1.  The  mathematical  complication  involved  in  a  changing 
mass  system,  while  possible  to  solve  on  the  computer,  would  require 
more  components  than  are  available. 

2.  The  mass  change  is  within  the  same  order  of  magnitude 
as  the  accuracy  of  the  computer,  this  there  would  be  no  observable 
gain  from  such  a  refinement. 

6.2  The  General  Differential  Equation  for  the  Internal  Combustion  Engine 0 
The  first  law  of  thermodynamics  may  be  written  in  differential 

form  as 

5  602.1 
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where  dQ  is   the  net  energy  entering  or  leaving  the  system  as  heat. 
dW  is   the  work  done  by  the   system. 
dE  is   the  change  in  internal  energy  of  the   system. 

By  taking  the  time  derivative  of  this  equation,   and  by  replaoing  dW  and 

dE  by  their  respective  equivalents  in  terms  of  pressure  and  volume,   the 

following  equation  is   obtained. 


6.2.2 


dt 

By  solving  for  the  highest  derivative  of  Pj 


and 


Now,  the  volume  is  given  by 


6.2.4 


6.2.5 


When  equations  6.2.4  and  6.2.6  are  substituted  into  equation 
6.2.2,  the  final  result  is 


*£#-*»*]■ 


6.2.6 


where  Pj;  is   the  engine  pressure. 

This  equation  may  be  solved  on  the  BEAC  in  the  following   steps. 

1.  Integrate     dPg  to  form  -PE  with  the  appropriate  initial 
oondition,   Pp. 

2.  Multiply  -PE  by  kx  to  form  -kPgx. 

3.  Add  k-1  dQ     to  form       Jill  d®   -    K PE  X\ 

AE     It  I  Ae     aft  J  ' 

4.  Divide  this  sum  by  x.      This  result  is,  by  equation  6.2.6, 

dPg,  which  is  then  fed  back  to  Step  1. 
W 
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The  above  steps  are  represented  schematically  in  Fig.  9.   ThiB  dia- 
gram also  shows  considerable  computer  detail  which  is  explained  in  the 
following  notes. 

1.  The  inputs  to  Amplifier  5  are  to  separate  the  net  dQ  into 

left 

"fuel"  and  "non-fuel"   terras  in  order  that  these  quantities  may  be  varied 

independently.      In  addition,   the  "fuel"  function  is  non-zero  for  only  a 

portion  of   the  cycle.      The   output  of  this  amplifier  is  dQ  net. 

7PE 

2.  The   relay  in  the  feedback  of  Amplifier  1  establishes  the 

initial  condition,  Pp,  whenever  x  is  greater  than  x., 

3.  The  combination  of  the  electronic  multiplier  and  Ampli- 
fier 2  is  pre-set  to  give  one-fiftieth  of  the  produot. 

4.  The  multiplier  in  the  feedback  of  Amplifier  4  has  the 
effect  of  dividing  by  x,  where -x  must  always  be  negative  for  stability. 
The  output  of  this  amplifier  is  pre-set  to  give   ten  times  the  quotient. 

5.  Amplifier  3  is  set  to  give  ten  times  the   sum  of   -  Jooi^^B^') 
anA  lh>  ("Tp"  ^?t  J    k*  order  to  compensate  for  the  constants  inherent  in 

the  multiplication  and  division. 

To  gain  a  better  insight  into   the  computer  solution  of  equation 
6.2.6  a  few  simple  qualitative  examples  will  be  discussed  in  the  follow* 
ing  three  sections. 
6.3  Computer  Simulation  of  a  Reversible  Adiabatic  Process. 

The  reversible  adiabatic  prooess  is  characterized  by  dQ     equal   to 
zero,   in  which  case  equation  6.2.6  reduces  to 

dP£  ,  KPEi 

It  may  be  easily  verified  that  Pgxk =  constant  is  the  solution  to  the 
above  equation,  which  defines  the  familiar  isentropio  process.   The 
oomputer  solution  to  equation  6.3.1  appears  in  Fig.  10  with  various 
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Figure  10o  Computer  Simulation  of  an  Isentropio  Proeese, 
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initial  conditions.  The  input  "x"  was  a  displaced  sinusoid  simulating  a 
piston  motion  at  about  one  cyole  per  seoond.  This  was  produced  from  the 
arrangement  shown  in  Fig.  11,  where  ^  is  an  auxiliary  variable  solving 
the  equation  f  +  UJ  j  -  O*  Also  shown  are  the  My"  and  H«N  variables 
whioh  will  be  used  in  the  forthcoming  discussions  of  the  compressor  and 
bounce  cylinders.  This  is  an  artificial  stroke  input  and  will  be  by- 
passed when  -the  entire  free  piston  analog  is  assembled.  *.  . 

The  solutions  depicted  in  Fig.  10  can  be  seen  to  be  excellent  for 
all  but  the  lowest  values  of  initial  oondition.   The  plotted  points  are 
the  numerical  values  of  the  indicated  equations,  and  the  solid  lines  are 
the  computer  solutions.   In  all  oases,  the  integration  was  stopped  when- 
ever "xn  was  greater  than  30  volts.   The  departure  from  the  correct 
solution  for  the  low  value  of  initial  oondition  is  due  to  the  fact  that 
the  inherent  electronic  errors  are  magnified  whenever  the  oomputer  is 
working  with  small  voltages.   This  oondition  must  be  avoided  if  at  all 
possible. 

It  is  to  be  pointed  out  that  the  multiplication  of  -Pg  by  x 
involves  two  quadrants.   That  is  to  say,  -Pg  is  always  negative,  while 
x  is  both  positive  and  negative.   Care  must  be  taken  to  insure  that 
this  multiplication  yields  the  same  result  in  both  quadrants.  A  slight 
maladjustment  in  this  multiplier  will  not  return  the  solution  to  the 

starting  point  since  dPg  in  one  direction  will  be  different  from  dPg  in 

dF  "cPT 

the  other. 

6.4  Computer  Simulation  of  a  Reversible  Polytropio  Process. 

The  reversible  polytropio  process  may  be  achieved  by  letting 


^(K-7l)PeX,  e^x 
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Upon  substitution,   equation  6,2.6  becomes 

4Pe   _        nPex 

—rrr     =      ~    6 '  6.4,2 

d£  X 

It  may  be  easily  verified  that  ?%xn   =  constant  is  the  solution  to  this 

equation.   It  should  be  pointed  out  that  while  mathematioally  dQ  is  as 

dt 
defined  in  equation  6.4.1,  the  same  result  may  be  achieved  on  the  oomputer 

by  letting  dQ  be  zero,  and  applying  the  coefficient  n  instead  of  k  after 

dt  2  "2 

Amplifier  2,  as  suggested  by  equation  6.4.2. 

Figure  12  illustrates  the  computer  solution  of  a  few  samples  of  this 

case,  again  with  excellent  agreement.  At  this  point,  it  is  important  to 

realize  that  there  is  no  net  loss  of  energy  from  the  system  over  a  cyole. 

The  energy  lost  during  compression  is  restored  in  the  same  manner  during 

the  expansion,  or 

This  may  be  easily  verified  on  the  oomputer  by  integrating  equation 
6.4.1.   Figure  13  shows  displacement,  pressure,  ^Q  ,  and  (/dcA^  as 

QaA  dt  ****' 

funotions  of  time.     As  expeoted,  ^dSa^tts  zero  at   the  end  of  a  oycle. 


6.5  Computer  Simulation  of  a  Non-poly tropic  Process. 

Consider  a  prooess  which  is  characterized  by  a  net  loss   (or  gain) 

of  energy.     This  may  be  accomplished  by  letting  dQ  be  some  funotion  other 

"oT 

than  that  defined  by  equation  6.4.1.     Suppose,   for  example,   that  jdQ 

dt 
were  a  small  negative  constant.      The  resulting  differential  equation  is 

not  easy  to  solve.     However,   it  oan  be  deduoed  that  a  oyclic  change  in 

volume  will  not  return  the  pressure   to  the  starting  point,   since     dPE 

w 

is  different  for  each  direotion  of  motion.  A  few  computer  solutions  of 
processes  of  this  type  are  shown  in  Fig.  14.  With  the  exception  of  the 
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Figurs  13o  Rersrsible  Pclytrcpis  Process  ss  a  Function  of  Tin*. 


31 


50 


? 

•H 


*  30 


20 


10 


0    u 


^ 


lc  20  3C  40 

Displacement,  rolts     (Arbitrary  Seal*) 


Figure  14o     Computer  Simulation  of  the  Nen-Polytropie 
Proeese0 


32 


isentropic  oase,  there  is  net  negative  work  done.   This  is,  of  course,  the 
work  equivalent  of  the  energy  lost  by  the  system  as  heat,  Figure  15 
illustrates  the  same  solutions  as  a  function  of  time.  Here  it  is  noted 
that  the  net  energy,  Q,  appears  as  a  ramp  function.   The  discontinuity  in 
the  pressure  curve  is  due  to  the  operation  of  the  initial  oondition  relay 
when  the  displacement  exceeded  30  volts. 

This  and  the  preceding  examples  serve  no  purpose  in  the  overall  aim 
of  this  investigation,  but  are  presented  in  order  that  the  reader  may 
beoome  familiar  with  the  rather  peculiar  sort  of  mathematical  manipula- 
tions used  in  setting  up  the  analog  computer.  In  addition,  they  may  serve 
as  useful  cheok  points  in  the  event  that  this  problem  is  set  up  on  a 
different  computer. 
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7C  Electronic  Analog  of  the  Standard  Diesel  Cyole. 

7.1  General 

Equation  6,2,6,  and  the  associated  computer  arrangement  illustrated 
in  Pig,  9  represents  the  most  general  oase  of  an  engine  oyole,  wherein 
any  form  of  heat  removal  or  addition  may  be  simulated  by  the  appropriate 
seleotion  of  functions  for  the  dQ  terms.   The  investigation  of  solutions 
of  this  type,  encompassing  various  methods  of  heat  removal  and  fuel 
injection,  is  a  very  broad  subject,  and  should  be  treated  as  a  separate 
investigation.   It  will  not  be  practical  to  do  this  within  the  wider  scope 
of  the  free  piston  problem.   In  order  that  the  results  from  this  investi- 
gation may  be  more  readily  analysed,  a  standard  diesel  cycle)  will  be 
employed  for  the  engine.  This  has  the  additional  practical  advantage 
that  excessive  computer  voltages  will  be  avoided0 

7.2  The  Standard  Diesel  Cyole, 

The  standard  diesel  cycle  consists  of  the  following  processes: 

1,  Polytropio  compression,  at  exponent  n,  of  the  engine  charge, 

2,  Constant  pressure  addition  of  energy  over  a  portion  of  the 
outstroke* 

3,  Polytropio  expansion  of  the  engine  charge  and  combustion 
products,  at  exponent  n' , 

4,  Constant  volume  heat  rejeotion  to  exhaust, 

5,  Constant  pressure  scavenging  process. 

Referring  to  equations  6,2,6  and  6,4,1,  these  processes  may  be  described 
mathematically  as  follows; 

1,  Compression  process. 

&t  —  - 
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and 


m-  *•$$»■■ 


and 


\  V     x/  7.2.2 

2.  Combustion  prooess. 

[A  o 

la't./p  /f-1  7.2.4  j, 

3.  Expansion  process, 
^^         =         -     *.'£*  7.2.5 


and 


%  -  *.(&»'■ 


7.2.6 


The  above  equations  may  be  readily  solved  on  the   computer  by  applying 
the  following  modifications  to  Fig.  9,  on  page  25. 

1.  Elimination  of  Amplifier  5. 

2.  Replacement  of  hk  by  |*i  when  the  velocity  is  negativej 
and  by  |n»  when  the  velooity  is  positive. 

3.  Disconnection  of  dPjg     from  Amplifier  1  during  the  pre- 

dt~ 

scribed  interval  of  oombusion,  thus  making  dPg  zero. 

dT 
Items  Z   and  3  above  may  be  accomplished  through  the  use  of  relays.   It 

should  be  noted,  however,  that  since  the  free  piston  engine  clearance  is 

variable,  the  constant  pressure  combustion  prooeBs  must  be  initiated  by 

a  ohange  in  sign  of  velocity  and  terminated  when  a  prescribed  amount  of 

energy  has  been  added  rather  than  at  particular  displacements.   The 

reason  for  this  is  that  a  change  in  clearance  should  not  affect  the 
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amount  of  fuel  energy  added. 

Figure  16  illustrates  the  computer  arrangement  to  simulate  the 
standard  diesel  cycle.  The  functions  of  the  additional  amplifiers  and 
relays  are  as  follows i 

1.  Relay  1,  actuated  by  a  ohange  in  sign  of  velocity*  applies 
the  coefficient  |n  when  velocity  is  negative,  and  the  coefficient  ^n'when 
the  velocity  is  positive. 

2.  Relay  2,   also  actuated  by  a  ohange  in  sign  of  velocity*  dis- 
connects dPg  from  Amplifier  1  when  the  veloolty  goes  positive!  and  connects 

dt~ 
the  circuit  so  that  dPg  may  be  reconnected  by  Relay  3  at  the  proper  time. 

dt~ 


3.  Amplifier  5  performs  the  integration  of  1^  n*PgX,  whenever 

10 
x  is  positive.   The  quantity  Jji'P-x.   may  be  rewritten  as 


This  integration,  therefore, yields  a  result  proportional  to  the  total 
instantaneous  fuel  energy.   This  is  compared  with  the  desired  value  of 
Qf  in  Amplifier  6  which  actuates  Relay  3  when   (/d0\  ^   equals  rn^HV. 

jUVn 

Figure  17  illustrates  a  typical   solution  using  this  arrangement.    It 
should  be  pointed  out  that  the  operation  of  this  analog  does  not  depend 
upon  a  pre-set  stroke  or   olearance. 
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8.  Energy  Relationships  in  the  Engine* 

In  the  operation  of  the  analog  of  the  standard  diesel  engine,  it  will 

be  necessary  to  know  the  net  energy  lost  as  heat  from  the  system;  and  the 

net  energy  added  during  the  combustion  process.   These  quantities  may  be 

determined  by  the  integration  of  equations  7.2.2,  7.2.4,  and  7.2.6  over 

the  appropriate  displacements.   These  arex 

1.   Energy  as  heat  transferred  during  compression 

-1 


8.1 


2,  Net  energy  as  heat  transferred  during  combustion 

3.  Energy  as  heat  transferred  during  expansion 


#*'  ~  ~rfw;  ^ Lw     "  ll  • 


8.3 


where  the   subscripts  1,2,3,  and  4  refer  to  the  points  illustrated  in 
Fig.   17.     It  can  be  seen  that  with  the  computer  furnishing  a  picture 
similar  to  Fig.   17,  all  of  the  required  pressures  and  volumes  will  be 
known,    thus  all  of  these  Q* s  may  be  determined.      The  sum  of  Qq  and  Qqi 
will  be  that  portion  of  Qrej  ohargeable  to  the  engine;   and  Qp  is  an 
idealized  rn^HV  product  which  does  not  include  heat  lost  during  combustion, 
nor  energy  lost  through  incomplete  combustion. 
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9,  Electronic  Analog  of  a  Reciprocating  Compressor. 

The  force  exerted  upon  the  compressor  cylinder  may  be  most  conven- 
iently represented  by  referring  to  a  typical  indicator  diagram  for  a 
compressor,  such  as  appears  in  Fig.  4(b),  page  9.   In  the  case  of  the 
free  piston  device,  there  is  no  fixed  geometrioal  point.  The  only  known 
quantities  are  the  pressure  limits. 

In  a  manner  similar  to  that  used  in  the  engine  equations,  the 
differential  equation  of  the  first  law  for  the  compressor  may  be  written 


as 


dt 


y  L 


Ac     dt  n/c^ 


9,1 


where  A.0  is  the  compressor  piston  area. 

Pc  is  the  compressor  pressure. 

y  is  the  compressor  displacement  from  the  head. 

dQ  is  the  rate  of  heat  transfer. 

dt" 
This  equation,  however,  is  valid  only  where  the  mass  is  constant,  not 

during  the  delivery  or  intake  processes.  For  these  processes,  it  can 

be  said  that 

P0  z  P(j»     for  delivery,  and 

9.2 
Po  =  pi#     ?or  intake. 

The  problem  of  the  changing  mass  system  of  the  compressor  may  be 

readily  solved  on  the  oomputer  with  essentially  the   same  arrangement  as 

that  used  for   the  engine,  but  with  the  use  of  diode  limiters  on  the 

amplifier  generating  P0.     Suoh  an  arrangement  is  shown  in  Fig.   18, 

which  solves  the  poly tropic  equation 


dt 


-TlcPc* 


9.3 
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It  oan  be  seen  from  Pig.  18  that  -P0  is  always  negative.  With  this 
being  the  case,  the  diodes  in  the  feedback  of  Amplifier  7  will  function 
as  follows t 

1.  With  y  negative  (inatroke),  dP_  will  be  positive,  and-P0, 

dt 
the  output  of  Amplifier  1  ,  will  tend  to  go  more  negative.   This  will 

continue  until-P0  overrides  the  bias  on  Diode  1  at  whioh  time  this  tuba 

oonduots,  thus  stopping  the  integration  with-Pc  held  constant  at-Pd. 

This  situation  will  continue  as  long  as  y  is  negative. 

2.  When  y  ohanges  sign,  dP  will  beoome  negative  and  the  out- 

dt 
put  of  this  amplifier  will  beoome  less  negative.  When  this  happens, 

the  initial  condition  will  be-Pd,  whioh  has  been  imposed  by  Diode  1. 

3.  The  absolute  value  of  P0  will  continue  to  deorease  until 
the  bias  on  Diode  2  overrides -Pc.   This  tube  will  then  conduct,  thus 
holding  the  output  at-P^.   This  will  continue  until  y  changes  sign  and 
the  integration  is  started  at  the  initial  condition -P^.  From  this  point 
on  the  oycle  will  repeat. 

Thus  it  oan  be  seen  that  by  the  use  of  two  diode  limiters,  equation 
7.1  may  be  solved  in  the  case  of  a  compressor  without  taking  into  con- 
sideration the  changing  mass  problem,  furthermore,  this  arrangement 
will  simulate  a  reciprooating  compressor  that  has  a  variable  stroke, 
whioh  is  precisely  the  situation  in  the  free  piston  maohine.   This  is 
beoause  the  integrations  oovering  the  compression  and  expansion  strokes 
are  initiated  by  a  change  in  sign  of  the  velooity,  not  at  a  particular 
position. 

The  analog  illustrated  in  Fig.  18  was  oonstruoted  using  the  assump- 
tion that  the  compression  and  expansion  prooeed  at  the  same  polytroplc 
exponent. 
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Figure   19  illustrates   two  solutions  of  the  compressor  problem, 
corresponding  to  two  different  stroke  lengths.     Here  it  is  again  noted 
that  a  pre-set  stroke  and  clearance   is  not  necessary. 
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10.     Energy  Relationships  in  the  Compressor. 

In  order   to  solve  equation  5.5.1  it  is  neoessary  to  know  the  amount 
of  energy  removed  as  heat  from  the  compressor  per  cycle.     Recall         that 
applying  the  coefficient  ntf  instead  of  k  is  equivalent  to  letting 

which  produces  a  reversible  polytropio  process  wherein  there  is  no  net 
loss  or  gain  of  energy  over  a  oyole.     This  will  be   true  only  in  the  case 
of  a  oonstant  mass  system.      The  compressor,   however,   compresses  one  mass 
and  expands  a  different  mass.      By  integrating  equation  8.1  over  the  com- 
pression and  expansion  strokes,    the  net  energy  transferred  as  heat  per 
oyole   is  . 

^  (K-l)(»c-l)  l[FJ        '  *J '    10.2 

where  m^  is  the  mass  of  air  delivered  by  the  compressor  per  cycle.   This 
equation  neglects  the  energy  transferred  during  the  delivery  and  intake 
processes.   It  can  be  seen  that  if  md  is  zero,  or  if  k  ■  no,  there  will 
be  no  heat  loss. 

*See  Appendix  III  for  derivation. 
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11.  Electronic  Analog  of  a  Oaa  Spring 

The  force  exerted  upon  the  bounoe  cylinder  is   that  of  a  gas  spring, 
i.e#)Pzk  r  constant.      In  this  case,  equation  6.3,1  may  be  rewritten  as 

dPb  /fP.Z 

-J £.    -      —      ,.L.  0      -      .  Ill 

dt    -  2        3  ' 

where  Pb  is  the  pressure  qxx   the  bounce  piston 
z  is  the  bounce  piston  displacement. 

On  the  basis  of  the  preoeding  analyses,  the  computer  arrangement 
may  be  immediately  depicted  as  that  shown  in  Fig.  20,  The  diode  in  the 
feedback  of  Amplifier  10  performs  the  function  of  an  atmospherio  relief 
by  not  permitting  Pb  to  fall  below  a  voltage  corresponding  to  one 
atmosphere,  thus  simulating  the  introduction  of    additional  air  into 
the  bounce  space,  A  typioal  solution  from  this  arrangement  is  illus- 
trated in  Pig,  21,  which  also  illustrates  the  functioning  of  the  at- 
mospheric relief. 

In  the  oase  of  the  bounce  cylinder  analog,  extreme  caution  must  be 
used  in  the  adjustment  of  the  electronic  multiplier  forming  the  Pz 
produot,  for  -the  reason  outlined  in  seotion  6.3.   This  situation  is 
aggravated  in  this  case  since  there  may  not  be  an  initial  condition 
imposed  every  cycle.   That  is  to  say,  the  solution  will  be  wfloating0M 
A  very  slight  maladjustment  in  this  multiplier  will  oause  the  solution 
to  drift  rapidly  away  from  the  correct  value.   This  malfunctioning  is 
illustrated  in  Fig.  22.   This  situation  must  be  avoided. 
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Figure  210  Cornuter  Simulation  of  af  Gee  Springc 
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12,  Coulomb  and  Viscous  Friction. 
12,1  The  Electronic  Analog. 

The  force  of  fr lotion  is  assumed  to  be  composed  of  the  two  familiar 
types:  Static  friotion,  a  constant  foroe  always  opposing  the  motion,  and 
visoous  friotion,  a  foroe  proportional  to  the  velocity,  again  opposing 
the  motion.  This  total  foroe  is  illustrated  graphically  in  Fig.  4(d). 
This  may  also  be  expressed  mathematically  as 

fy   -  C  "  nri      (  i  <  c) 


Ff   =  -C-*r'i      {i>o)  > 

where  C  is  the  Coulomb  friction  force 

v  is  the  viscous  friction  coefficient. 

The  analog  to  produce  a  voltage  corresponding  to  this  foroe  is 
illustrated  below  in  Fig.  23. 


12.1 


-www- 
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Figure  23.  Electronic  Analog  of  Coulomb  and  Viscous  Friction 
The  amplifier  has  no  feedback  except  for  the  two  diodes  and  their 
series  resistances.  Whenever  the  input  is  different  from  zero,  even  by 
the  smallest  amount,  the  amplifier  output  will  tend  immediately  towards  ■ 
saturation.   This  tendency  is  checked  at  the  voltage  "C  by  one  or  the 
other  of  the  diodes,  which  when  conducting,  places  the  feedback  resist- 
ance in  the  circuit.  Figure  24  illustrates  the  computer  analog  of  the 
foroe  of  friotion  both  as  a  function  of  velocity  and  of  displacement. 
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12,2  Energy  Dissipated  by  PriotionQ 

The  last  remaining  item  as  yet  unknown  in  equation  5.5,1  is  that 
portion  of  Qrej  chargeable  to  friotion. 
The  force  of  Friction  is  given  by 

'p  ~   C  ~  VX   j  when  x  is  negative, 

and  12.2.1 

h      ~     ""  C  ""  V~X-  j  when  x  is  positive. 

In  order  that  this  force  may  be  integrated  directly,  it  shall  be  assumed 
that  the  free  piston  has  the  approximate  motion 

X   »  Ce  t    ^-/c**,  J7rit  ^  12.2.2 

where  oe  is  the  engine  clearance 
f  is  the  frequency 
s  is  the  stroke 
The  work  done  by  friotion  is 


'fm       §    F+cLf., 


or  JL 


wF*  S  f^)^- 


12.2.3 
12.2.4 


O 
Upon  substituting  and  integrating  over  a  cycle,  the  net  energy  lost  through 
friotion  work  is  given  by 

Vl/c   -        3CS      +       12L£^  12.2.5 

r  a. 
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13.  Computer  Operated  Relays. 

As  outlined  in  Pig.  16,  the  use  of  five  relays  is  required  in  the 
simulation  of  a  diesel  engine.  For  convenience,  these  are: 

10  A  relay  to  start  and  stop  integration  whenever  x  is  greater 
than  xp,  the  port  position.  This  has  the  effect  of  putting  the  "tail"  on 
the  engine  picture,  thus  simulating  the  scavenging  prooess. 

2.  Three  relays  necessary   to  make  dPg  zero  for  a  portion  of 

dt" 

the  outstroke,  thus  simulating  the  oonstant  pressure  addition  of  fuel. 

3.  One  relay  to  ohange  exponent  from  compression  to  expansion. 
The  above  relays,  of  course,  must  be  actuated  by  quantities  available  in 
the  computer.   In  all  oases  the  relays  were  driven  by  the  output  of  an 
operational  amplifier  by  use  of  the  scheme  illustrated  below  in  Fig.  25. 


et   #-A4MV-i 
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Figure  25.  Use  of  Operational  Amplifier  to  Aotuate  Relays. 
As  can  be  seen  above,  the  amplifier  has  no  feedbaok  except  diodes  to 
limit  the  output  to  a  voltage,  £jSuff icient  to  actuate  the  relay.  The 
action  of  the  above  arrangement  is  to  yield  an  immediate  step  voltage  to 
the  relay  whenever  the  algebraic  sum  of  the  inputs  passes  through  zero0 
Figure  26  illustrates  the  output  of  such  an  amplifier  actuated  by  a 
sine  wave.  It  can  be  seen  that  this  output  is  very  nearly  a  square 
wave0 
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14,  Electronic  Analog  of  the  Free  Piston  Generator, 

The  preceding  sections  may  be  summarized  by  referring  to  an  overall 

computer  arrangement  for  the  solution  of  the  piston  motion.   This  is 

illustrated  in  Fig.  27.   This  diagram  is  basically  the  same  as  Fig.  7, 

except  that  the  details  of  Figs.  16,  18,  20,  23,  and  25  are  included. 

The  numbers  on  the  amplifiers  correspond  to  those  used  in  the  above 
mentioned  figures. 
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16.  Time  and  Magnitude  Sealing, 

16.1  General. 

All  of  the  preceding  seotions  have  been  rather  qualitative  in  nature 
in  order  that  the  principles  of  the  analog  approach  to  the  problem  may  be 
presented  more  clearly.  In  order  to  obtain  a  numerical  result,  however, 
it  is  necessary  to  consider  magnitude  and  time  soaling  within  the  computer. 

15.2  Soaling  the  Pressure  Analogs. 

As  outlined  in  section  6.3,  the  various  pressure  analogs  were  con- 
structed by  using  as  a  driving  force  a  separately  generated  displaced  sine 
wave,  operating  at  a  fixed  frequency,  based  on  real  time.  Such  will  not 
be  the  case  when  each  of  these  analogs  becomes  a  part  of  the  overall 
hook-up  as  illustrated  in  Fig.  27.   In  this  case,  the  integrations  form- 
ing displacement  and  velocity  will  seek  some  natural  frequency,  based  on 
computer  time,  not  real  time.   To  this  end,  a  time  soaling  factor  of  100 
will  be  used.   Thus,  a  real  machine  operating  at  600  cycles  per  minute 
will  appear  on  the  computer  as  operating  at  6  cycles  per  minute.  The 
only  criterion  necessary  in  making  the  pressure  integrations  compatible 
time-wise  with  the  displacement  integrations  is  that  a  common  time  scale 
be  used.  When  this  is  done,  all  integrations  will  proceed  at  the  same 
rate. 

The- magnitude  scaling  is  entirely  arbitrary,  depending  upon  the 
maximum  value  of  the  variable,  and  its  relationship  to  the  50  volt  limi- 
tation of  the  BEAC  amplifiers.   The  only  criterion  here  is  that  the  same 
soale  factor  must  be  applied  to  the  initial  conditions. 

There  is  a  definite  relationship  between  the  various  soale  factors 
involved  within  each  computing  loop  simulating  the  engine,  compressor, 
and  bounce  pressures.  The  compressor,  for  example, solves  the  equation 
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P  =  ftfk  =    -  nc  Pc   j 
**         dt  y  "  '*  15.2.1 

or,  in  terms  of  computer  voltages 

where  Q^  is  the  magnitude  scale,  units  of  "iM  per  volt,  and  the  barred 
quantities  represent  the  oomputer  voltages  corresponding  to  the  real 
quantity.  By  comparing  the  above  equations,  it  can  be  seen  that  C^o 
may  not  be  arbitrarily  assigned,  but  is  given  by 

<Xp      =    J^&J^L  .  15.2.3 

Similar  expressions  exist  in  the  engine  and  bounce  analogs.  These  are 


Oi 


1      *  °^°$ 


Pe  *_ ^  15#2o4 


and 


OCRt     =   ll *-  ,  15.2.5 


ra 


<** 


As  outlined  in  section  6.2,  there  is  no  liberty  selecting  the 
magnitude  soales  involved  in  the  multiplications  and  divisions.  This 
scaling  has  already  been  performed  for  all  amplifiers  except  those  in- 
tegrating to  form  the  various  pressures. 

In  the  case  of  the  engine,  for  example,  the  integrating  amplifier 
may  be  treated  as  follows  t 
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Xp 


p    a, 


-P. 


In  this  amplifier, 


p    3 


15.2.6 


where  CXt  is  the  time  scaling  factor,  computer  time. 

real  time 

When  the  above  constants  are  associated  with  the  circuit  parameters  as 

outlined  in  Appendix H,  the  following  equation  is  formed 


OTp 


Oi0    <X 


Pp  ^t 


15.2.7 


#/ 


When  equation  15.2.3  is  substituted. 


a* 


c*  ^ 


/?,  C 


*i 


c*-,  c* 


15.2.8 


*  ^Z- 


By  similar  arguments,  the  circuit  parameters  involved  in  the  com- 
pressor and  bounce  pressure  integrations  are   given  byi 


and 


&7    .     m 


2Lh. 


o^  c** 


15.2.9 


& 


ID 


Rio  C 


to  ^io 


15.2.10 
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The  subscripts  1,  7,  and  10  used  in  the  above  equations  refer  to  the 
number  of  the  amplifier  as  given  in  Fig.  27. 

Since  x,  y,  and  z  refer  to  the  same  motion,  it  must  be  true  that 


o(x  =r    cyM    -    <* 


?  > 


15.2,11 


and 


<**  -  «5  -  «i 


15.2.12 


When  this  substitution  is  made,  there  finally  results 

a,  Sly  a.,0 


Ft,  C*i  fty  Cz 


7^7 


n,n  C. 


10  *-$io 


15.2.13 


This  relationship  specifies  the  definite  time-magnitude  relationship 
between  the  pressure  integrations  and  the  displacement  integration. 
15 o 3     Scaling  the  Displacement  Integrationso 

The     displacement  integrations  may  best  be  analysed  by  redrawing  a 
portion  of  Fig.  27: 


-ft 


-iw 


+S±t2L 


0^/13 


0dm 


The  following  may  be  written  for  Amplifier  H. 


o$X    = 


-[-  /w/*  -  A^£  *  /^/j>  +  *v£j;  15- 


3ol 
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From  which 


&i+l]\. 


K 


Jlfl      _  /if  vt-Pe 


As  OCt 


HI 


M  ex 


15.3.2 


15.3.3 


and 


A/*3    ^W3  nfe   ** 


^ 


fl, 


f3 


H  <*, 


*/^ 


C* 


/* 


/W 


//      <*v 


15.3.4 


16.3.5 


The  following  may  be  written  for  Amplifier  15 


<***  -   - 


sauLdt 


or 


<2 


/£ 


C*y 


R\S  CjA5 


ctf*   <*■* 


15.3.6 


15.3.7 


And  for  Amplifier  l6j 


cx±  J 


15.3.8 


or 


<2./6 


#/6  4wfc 


o<*  <** 


15.3.9 


By  comparing  this  equation  with  equation  15.2.13.  it  can  be  seen  that 
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L&_  . 


otx&t 


a7 


A  c«         /?-  C 


/* 


7  ^-&7 


•    15.3.10 


This  equation  means  that  the  RC  time  constants  in  the  x,  PE#  ^  ,  and  F^ 
integrations  must  be  the  same. 
*  15.4  Scaling  the  Friction  Analog. 


Again  referring  to  a  portion  of  Fig.  27 i 


T^y 


X       3L 


13 


5 


The  equation  for  this  amplifier  is 


<yFiF&  -;abC  -    vo^x  J 


15.4.1 


from  ishich 


an  ft 


^13       _         Ofr  V 


/?, 


15.4.2 


No  further   scaling  is  required  in  the  other  amplifiers  of  Fig.   27. 
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16.  The  Pescara  Machine  on  the  Analog  Computer. 
16.1  Results. 

Reference  1  contains  the  data  and  a  digital  analysis  of  a  free  piston 
gas  generator  manufactured  by  the  Pescara  organization  in  France.  Appen- 
dix IV  contains  the  applicable  original  data  and  some  of  the  results  given 
by  London.   It  is  interesting  to  note  that  the  analog  approaoh  to  the 
problem  is  nearly  the  reverse  from  the  digital  approach.  The  digital 
approaoh  uses  as  a  starting  point  an  assumed  stroke  and  clearances,  and 
then  prooeeds  to  caloulate  some  of  the  geometry,  frequency,  power  output, 
and  other  quantities  of  interest.   The  analog  computer  approach  uses  as 
a  starting  point  a  particular  geometry,  and  yields  as  a  solution  stroke, 
clearance,  and  frequency.  Unfortunately,  the  data  required  by  the  com- 
puter were  not  entirely  complete  in  the  original  Pescara  paper. 
Conspicuously  absent  were  J 

1.  The  piston  mass,  M. 

2.  The  absolute  bounce  clearance,  b. 

3.  The  initial  bounce  pressure;PD. 

All  other  items  were  either  given  or  were  determined  by  simple  calcula- 
tion. 

Of  oourse,  the  computer  cannot  produce  a  solution  in  the  absence  of 
the  above  items.   London,  however,  had  calculated  these  items  in  his 
analysis,  and  his  results  being  the  best  available  data,  were  used  for 
these  three  parameters.   It  should  be  pointed  out,  however,  that  this 
digital  analysis  is  not  quite  compatible  with  the  computer  solution  in 
thatt 

1.  The  digital  analysis  used  a  modified  Otto  engine  cycle, 
while  the  computer  uses  a  diesel  cycle. 

2.  The  digital  method  oharges  friction  individually  to  the 
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engine  and  compressor  in  the  form  of  meohanioal  efficiencies,  while  such 
is  charged  to  the  maohine  as  a  whole  in  the  computer  approaoh. 

3.  The  digital  method  uses  "lumped"  quantities  such  as  net 
engine  work,  net  compressor  work,  and  net  friction  work.  Such  is  not 
possible  on  the  computer  since  the  computer  must  consider  differential 
items  in  solving  the  equation  of  motion. 

With  these  exceptions  noted,  the  following  compares  the  oomputer 
solution*  of  the  Pesoara  maohine  with  the  combined  original  data  and 
London' 8  analysis:* 

1.  Geometrical  Results 
Item 

Frequency,  cycle s/min. 
Stroke,  inches 

Engine  effective  stroke,  inches 
Engine  clearance,  inches 
Compressor  clearance,  inches 
Bounce  Clearance 

2.  Mass  Results 


Compressor  air  delivered 
lbs/side/cyole 


Analog 

Pes car a 

600 

613 

17,5 

17.6 

9.74 

9.70 

1.10 

1.14 

2.33 

2.38 

2.54 

2.50 

Ratio,  eagine/comp,  air 
Air-fuel  ratio 
Gas  flow  rate,  lbs./hr. 
Fuel  flow  rate,  lbs./hr. 


0.500 

0.580 

0.370 

0.516 

22.7 

32.0 

36,600 

29,000 

589 

460 

Refer  to  Appendix  17  for  scaling  of  Pescara  Machine 
"Results  and  Calculations  appear  in  Appendix  7. 
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Analog 

Pesoara 

1550°R 

1405  °R 

1200 

1138 

3.   Thermodynamic  Results. 
Item 

Delivery  Temperature 
Air  Horsepower 
The  net  energy  balance  was  also  calculated.  This  showed  a  discrepancy 

of  6.6^  between  the  engine  work  and  the  compressor  plus  frictional  work0 
Figures  28,  29,  and  30  are  the  computer  solutions  of  the  engine 

pressure,  compressor  pressure,  and  bounce  pressure,  respectively,  as 

functions  of  their  respective  displacements. 

16.2  Discussion  of  Results. 

As  shown  in  the  preceding  seotion,  the  geometrical  results  were 

in  excellent  agreement  with  the  actual  machine  data,   'The  rest  of  the 

results,  however,  were  close,  but  not  as  good  as  might  be  expected. 

This  is  attributable  to  two  causes: 

1.  The  analog  computer  was  actually  solving  a  different  kind 
of  problan  than  the  Pescara  machine,  even  though  the  geometry  was  essen- 
tially the  same.   This  is  because: 

a.  A  diesel  cycle  was  used  instead  of  a  real  engine  cycle. 

b.  Different  thermodynamic  constants  were  used. 

c.  The  frictional  constants  were  assumed  at  much  too  low 
a  value. 

2.  Computer  error.   In  this  regard  it  is  to  be  noted  that 
extreme  difficulty  was  encountered  in  making  the  bounce  analog  stable. 
This  phenoraonon  resulted  from  the  error  in  forming  the  Pz  product  by  an 
electronic  device.  Since  z  changes  sign  it  is  imperative  that* 

pz  ■  -p(-z)f  so  that  the  solution  will  reproduce  itself  regardless  of 
direction  of  motion.   The  electronic  multiplier  is  not  capable  of  the 
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Figure  2S0     Computer  Sisruletion  of  the  Engine  of  the 
Pe scare  M?chine0 
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Figure  29o  Computer  Simulation  ol  the  Compressor 
of  the  Pe  see  ret  M*tehi»c0 
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Figure  30o  Computer  Simulation  of  the  Bounce 
of  the  Peseara  Machine 0 
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accuracy  that  is  required  in  the  case  of  the  bounce  piston.   The  engine 
and  compressor  analogs  are  not  sensitive  to  this  effect  in  that  an  initial 
condition  is  imposed  each  cycle  which  prohibits  the  solution  from  drifting 
away.   Thus,  the  analog  of  the  bounce  piston  exhibited  a  phenomenon 
similar  to,  but  not  as  pronounced  as,  that  shovnn  in  Fig.  23,  page  51. 
The  maximum  drift  was  of  the  order  of  3  volts.   The  scaling,  however 
was  such  that  a  3  volt  departure  represented  a  force  error  of  about 
59,000  pounds,  which  i3  sufficient  to  disrupt  the  entire  analog,  causing 

the  solution  to  become  unstable  after  about  three  cycles.  Thus  the 

i 

solution  as  outlined  in  Appendix  V  approaches,   but  may  not  be,    the  steady 
state.      This  effect  may  account  for  most  of  the  discrepancies. 
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17.   Conclusions 

The  free  piston  engine  problem  may  be  solved  on  an  eleotronio 
analog  oomputer,  if  such  a  computer  is  capable  of  generating  a  stable 
solution  of  the  gas  spring.   The  computer  can  produce  a  solution  to  the 
equation  of  the  first  law  of  thermodynamics  in  the  simulation  of  an  in- 
ternal combustion  engine  and  of  a  reciprocating  compressor.  An  analog 
computer  is  not  capable  of  solving  the  gas  spring  problem  to  the  required 
aocuracy  in  this  manner  due  to  small  inaccuracies  in  electronic  multipli- 
cation. 

Two  alternative  ways  are  suggested  to  produce  an  accurate  analog 
of  the  gas  spring.   These  arei 

1,  Use  of  a  function  generator  to  yield  a  funotion  which  may 
be  used  to  give  Pz^  -   constant.  Such  a  device  would  keep  positive  control 
of  the  bounce  force  at  all  times, 

2.  Use  of  one-quadrant  multiplication.   It  has  been  pointed 
out  that  the  primary  difficulty  was  that  electronic  multiplication  could 
not  give  exactly  Pz  s  -P(-z).   This  difficulty  may  be  avoided  by  using 

z  having  always  the  same  sign.   This  may  be  accomplished  through  the  use 
of  velocity  relays  and  a  sign  changer  to  switch  both  the  input  and  out- 
put of  the  multiplication  such  that  the  multiplier  operates  only  in  one 
quadrant. 

It  should  also  be  pointed  out  that  the  solution  of  this  problem 
required  a  good  sized  computer,  in  that  a  total  of  26  amplifiers,  eight 
diode  limiters,  eight  voltage  sources,  and  six  electronic  multipliers 
were  required,  in  addition  to  the  necessary  read-out  equipment.   This 
is  considered  the  minimum  computer  size  to  handle  suoh  a  problem. 
Further  refinements,  or  the  addition  of  another  compressor  analog  for 
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the  solution  of  the  free  piston  compressor  would,  of  oourse,  inorease 

the  computer  requirements. 
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APPENDIX  I 
LIST  OF  ASSUMPTIONS 


1.   The  scavenge  air  is  furnished  at  constant  pressure  and 
temperature  o 

20   The  delivery  to  the  turbine  occurs  at  constant  pressure 
and  tempera ture0 

3o   The  compressor  intake  is  70°Fj  one  atmosphere0 

4o   The  operating  fluid  obeys  the  ideal  gas  laws* 

5o   The  pressure  drop  across  each  valve  and  engine  port  is 
five  per  cento 

60   The  geometry  is  symmetrical  with  respect  the  centerline0 

7o   The  compressor  compresses  and  expands  in  accordance  with 
a  polytropic  exponent  of  1.30. 

80   The  sensible  enthalpy  of  the  fuel  may  be  neglected, 

9o   The  engine  operates  in  accordance  with  an  air  standard 
diesel  cycle. 

10o   The  energy  transferred  as  heat  during  the  compressor  intake 
and  delivery  is  neglected. 

llo   The  frictional  force  consist*  of  coulomb  and  viscous 
friction 0 

120   The  friction  work  is  determined  in  advance  by  assuming 
that  the  stroke  is  sinesoidalo 


Other  numerical  assumptions  regarding  the  specific  problem 
of  the  Pescara  machine  are  found  in  Appendix  IV. 
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APPENDIX  II 
Analog  Computer  Symbols 


a) 


Circuit 


Symbol 


Function 


1«  Potentiometer: 

e, 

e9 


e, 


e,  a  <a  e, 


2*  Summer: 


ft 

r-A/VWV-r 


1/ 


S— +22- 


3«  Summing  Integrator: 


H 


^ 


<5 


4.  Relay  Driver: 


e« 


,£, 


5> 


> 


^ 


■-^♦^f* 


*   <3j% 


lE"e* 


* 


*.~m**Jtd 


0 


A, 


A"3' 


^     a*    fc 


% 


*Uft 


&*£*-_ 


2_ 


3 


L3 


At- 6   wAe-n  Jill  +  af€i  <  ft 
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APPENDIX  II   (Continued) 

Circuit  Symbol 

5o  Beeimg  Electronic  Multiplier x 
10 


Meaning 


e, — | 


BbM   ■* 


ft  r+lfe  Q  ** 


*>3  '     3D 


Notes     Direct  multiplier  output  is  approximately    -      5  fi*  _ 
60  Donner  Electronic  Multiplier: 


DF/y 


<??=  ~  /0#  ^/^a 


70  Diriaion  using  Boeing  Multiplier: 


a, 


r  Bf« 


ea 


r-fif/f 


e* 


e,  a, 


e*  ft       ea 


80  Division  Using  Donner  Multiplier: 


DFAf 


eA 


3 


a.   n 


^ 


ea  *  -/*£> 


Other  Special  Circuits  are  explained  in  textc 
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APPENDIX   III 
Ifeat  Removal  from  the  Compressor 
The  compressor  cycle   is  given  by  the  figure  below: 

e 


For  the  processes  1-2  and  3-4,  the  rate  of  heat  transfer  is  given  by 


It  shall  be  assumed  that  the  delivery  and  intake  processes  involve 
no  energy  transfer  as  heate 
For  process  1-2 


p-  /:(*) 


J 


or 


Upon  integrating  between  the  volume  limits  1  and  Z, 
/O   -   A     ^'  ^c  P   ""c  I 


l  y 


*»«-■ 


&*•"'/> 


Now 


and 


4y,--v;3 


«*U-  (*) 


««-*' 


py,--  tk,  a  % . 


Therefore 


Q, 


K-r\t 
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Rtk 


>{**-  V). 


By  a  similar  process,  Q34  is  given  by 


Qi 


s-.,  R»,  IT.  -  71 


The  net  energy  transfer  is 

5,  -    £«  +  a 


f 


K~*< 


now 


i*        ^3  Y 


^  -*,  --   -  7hdL. 


[K-l)l»e-l 


(T-T)(?ns-^)FL3 


Therefore 


&  - 


_   /f-n 


{K-i)(-r>i- 


ij^^[f  -1 


1  ^V^^ 


1 
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APPENDIX  17 

PROGRAMMING 

THE  ANALOG  COMPUTER  FOR  THE  PESCARA  MACHINE 

I.  Data. 

Listed  below  is 

the 

combined  data  from  the  original  data  and  the 

results  given  by  London t 

(1) 

A.  Geometry* 

A« 

0o979  ft2   (*) 

Ac 

5.73  ft2   (*) 

Ab 

6.77  ft2   (*) 

M 

2080  lb    (#) 

XP 

10.84"      (f) 

0 

1.24"      (t) 

b 

21.14"      (t) 

B.  Operating  Conditions t 

P*  4.85      (*) 

P-b  64.5psia    (*) 


* 

293  psia    (£ 

c. 

Assumed  Quantities! 

C 

50  lb 

V 

0.1  lb  sec/ft 

n 

1.30 

w  nt 

1.40 

n0 

1.30 

♦Original  Data. 

f Calculated  from  original  data. 

+  From  London's  solution. 
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II,  Assignment  of  Scale  Paotors. 

A,  Pressure  Soaling. 

1,  The  Engine, 

Pp  z  64.5  psia 

The  compression  ratio  is  estimated  as  10 
Therefore   Peak  pressure    ft*  1500  psia 
Let  ot^   s  30  psia/volt  a  4320psf a/volt 
Then  Pp     a  2,15  volts, 

2,  The  Compressor, 
P±  a  14,0  psia 

Pd  ■  14,0  x  4.85  •  67,9  psia 

Lot  o^/»  a  10  psia/volt  =  1440  psfa/volt 

P^  -  1,40  volt 8 

Pd  s  6,79  volts 

3,  The  Bounce. 

Let  otp    =    20  psia/volt  2  2880  psfa/volt 
Pg  •  14,65  volts 

B,  Displacement  Scaling. 

Let  o^    «  |n  per  volt  ■  l/24  ft  per  volt 
Then^Xp  •  21.68  volts 

o"    =  2,46  volts 

b"    s  42,28  volts 

C,  Velocity  Scaling, 
It  is  knoira  that 

x«sl  -|-    9  cos  62. 8t    inches 

*max  m  50       ftA«0 

Let  06    «  2  ft/seo/volt 
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D.  Acceleration  Soaling. 

xmax  *K  3000  ft/seo2 

let  Otj   »  100  ft/seo2/Wt 

E.  Time  Soaling. 

let  °<£  a  100  seconds  computer  time/l  seoond   real  time 

F.  Friction  Foroe  Soaling. 

let  <XF  -  20  lbs/volt 

then  C  =  2.5  volt  s 
0.     Summary. 

°^PF  a  4320  psf a/volt 

0£pc  s  1440  psf  a/volt 

°^  =  2880  psfa/volt 

O^x  =  1/24  ft/volt 

cyx  =  2  ft/seo/volt 

<^y  =  100  ft/seo2/volt 

&FW  =  20  lb/volt 

°^t       s     100  seo/sec 
III.   Determination  of  Ciroult  Quantities 

A.     Pressure  and  Displacement   Integrators   (Amplifiers  1,   7,   10,   and 
16.) 

From  equation  15,3. 10 j 


R,C^,         )%C^y  Rio  C$io  R/&  C&IC  oCx  CXt 


=    0.480. 


Let  R's  -  1  meg 

C»S    *    1    M? 

a's  •  0.480 
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B.  The  velocity  Integration  (Amplifier  15) 
From  equation  15.3,7 

<%/$ ._    =      — ^&~-   -       0.500 

Let 

C.  Summation  of  Foroes   (Amplifier  14) 
From  equations  15,3.2  to  15.3.5.' 


<£,V/    R 


$1+ 


ft 


If  I 


S-H2    R 


*IH 


Q.C&4- 

Me*; 

±279 

Ag^ 

3.0  SL 

a'^  Rk'r       -  o^>       ^  0.0031 

Lett 

/f/v,?    =    *5  Key  3  i$J2/n~    <$-/<+#  -  Oo  39 

ft,+3     =     jTvejjTfo**.     £,f3    -  0.3O2 

$w    -      /0  ^   2x*~    £•/*?    -  0,03/0. 
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D.  Formation  of  Frictional  Force  (Amplifier  13) 
From  equation  15,4.2 


3-i  3  K$n 

A 13 

<** 

O.CL 

&t: 

R*3    - 

OJ  Tn^y 

tf«  - 

10     ?*efy  ,    ~t*M^ 

al3  =  100 
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APPENDIX  V 
CALCULATION  OF  RESULTS  FROM  THE  ANALOG  COMPUTER  SOLUTION 


1*  Frequency, 

The  measured  frequency  was  six  cycles  per  computer 
minute,  corresponding  to  a  real  frequency  of  600  cycles 
per  minute. 

2.  Stroke • and  clearances. 

Stroke  from  figure  28  -  17.5  inches. 
Engine  effective  stroke  -  9.74  inches. 
Engine  clearance  -  1.10  inches. 
Compressor  clearance  -  2.35  inches. 
Bounce  clearance  -  2.54  inches. 

3.  Compressor  delivery. 

From  equation  5.5.2:  m<j  »   PA  AQ  Yj/  RT     . 

1    I      - 

Y,  *    29.9  volts  «=  14.95  in. 

P*  «    14.0  psia 

A*  «      5.73  ft2 

Tq  *   530°  R 


*  ■  {u<%immi?' 


95)     =   0.500  lb 


4.  Engine  charge  mass. 

From  ecuation  5.4.5 i 


me/Vp  *  0.0113  Pt0#7 
■  0.209  lb/ft3 

V  » (10.84) (0.979)  «   0.884  ft3 

P        12 


83 


4.  (Continued) 


in  s  (0,209)  (0.884)  =  0.185  lb, 

9 


5.  Ratio  of  compressor  to  engine  air< 

me/md  ss  0.185/0.500  «  0.370. 

6.  Mass  of  fuel. 


From  equation  8.2  :   mfHV»  k  P2  (V3*-V2) 


F2  » 51  volts  «  1530  psia 
SV  »  3.6  volts  »  1.80  in. 

nLpHV  *  1.4  (1530) (144) (1.8)  =115,500  ft-lb 
^     ff^  (I2~) 

Assuming  a  LHV  of  18,200  BTU/lb; 

hlp«   115.500     ~  0.00818  lb. 
f   (13.2$0)(^8) 

7.  Air-fuel  ratio. 

a-f  ratior  0.185/0.00818  w  22.7 

8.  Fuel  rite. 

Fuel  rate*(0.00818)(2)(600)(60)  -  589  lb/hr 

9.  Gas  rate. 

Total  mass  per  stroke  per  side  «  0.508  lb. 

Gas  rate*  (0.508)  (2)  (600)  (60)  *  36,600  lb/hr 
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10.  Friction  work. 

From  equation  12.2.5;  Wf  -  2Gb  ^-i^rfB2, 


Wf=-  (2)(50)(17^5)  ,  (7T*)(0ol)(17^5)2   =:  153  ft-BD 


11»  Heat  lost  "by  compressor. 


From  equation  10.2,   Q.   *  J^JH-yQ  fill  f  f  p*)^t       A~] 


=  5170  ft-lb. 

12 .Heat  lost  from  engine. 

From  equation  8.1,       Q    »     fr-*J  % V'        /KT  _    4 

*"        fc&-0   [  I  If/ 

Q©   -     (0.1) (0*979) (64.5) (144) (10.84)     /"  n  „ 

6  (0\4)f6.3Hi2)  f-  (8.5)0*3  -  1 

«      8,000  ft  lb 

13.    Temperature  of  exit  gasses. 
From  equation  51  5.5.1: 

mdh0+  mjHV  =■  (mf-t-md)^  +  Qrej 

Bid h    =  (0.500)  (778)  (126.7)        49,200  ft-lb 
mfHV'-  115,500  ft-lb 

^rej  =  13,300  ft-lb 

from  which 

"» -  tgfaffl.w>  "  362Bro/lb 

From  the  tables? 
Tt  -  1550OR# 
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14.  Air  Horsepower. 

From  the  gas  tables,  an  isentropic  expansion 
from  the  state  given  in  item  13  yields  an 
ideal  work  output  of  139  BTU/lb. 


AHP  _ (159) (778) (0.508) (2) (600) (60)    -   1200 

33,0000 


15,  flork  balance. 

i 

Area  under  the  engine  curve  "340  volts  • 
We=-( 340)  (4320)  (0.979)  (1/24)  -=  60, 000  ft-lbs 
krea  under  the  compressor  curve  -=  160  volts2. 
Wc= (160) (1440) (5.73) (1/24)  =     55,200  ft -lbs. 
Wf  =  153  ft-lbs. 

Discrepancy:  4,650  ft-lbs 


! 
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APPENDIX  VI 
THE  BOEING  ELECTRONIC  ANALOG  COMPUTER 
The  Boeing  Electronic  Analog  Computer  (BEAC)  consists  of  an  equip- 
ment rack,  in  which  is  mounted  a  number  of  operational  amplifiers, 
coefficient  potentiometers,  iimiters,  multipliers,  and  voltage  sources. 
The  rttok  also  contains  a  common  power  supply  for  the  various  items.  For 
this  investigation,  two  such  raoks  were  used.  A  photograph  of  this  equip- 
ment is  shown  in  Fig.  31.  Also  illustrated  is  the  x-y  plotting  device, 
and  the  time  reoorder.   The  hook-up  to  solve  a  particular  problem  may  be 

made  in  two  wayst 

I  i 
1.  Direotly  on  the  front  of  the  computer,  or 

20  On  a  pre-patch  board,  which  is  then  mounted  on  the  front  of 

the  machine.  All  the  necessary  connections  are  led  in  beneath  the  pre- 

patoh  panel. 

The  toeing  operational  amplifier  is  a  high-quality  DC  amplifier 

having  an  open  loop  gain  of  approximately  50,000.   This  amplifier  will 

operate  with  stability  with  any  feedback  ratio  up  to  unity.   It  .may  be 

used  for  summation,  integration,  differentiation, or  more  oomplex 

operations.   The  amplifier  remains  linear  as  long  as  the  output  is  less 

than  50  volts,  A  glow  tube  circuit  shorts  the  output  at  about  51  volts, 

thus  this  value  is  an  absolute  operating  limit0 

The  general  operation  of  the  oomputer  is  quite  stable.  This 

i 

:| 

stability  may  be  enhanced  further  through  the  use  of  an  automatic  bal- 
ancing system,  whioh  restriots  the  amplifier  drift  to  less  than  a  few 
millivolts  per  month. 

Detailed  descriptions  and  operating  procedures  may  be  found  in  the 
manufacturer's   data,v    '< 
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